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We study the random walk problem on a deterministic scale-lree network, in the presence of a 
set of static, identical targets; due to the strong inhomogeneity of the underlying structure the 
mean first-passage time (MFPT), meant as a measure of transport efficiency, is expected to depend 
sensitively on the position of targets. We consider several spatial arrangements for targets and we 
calculate, mainly rigorously, the related MFPT, where the average is taken over all possible starting 
points and over all possible paths. For all the cases studied, the MFPT asymptotically scales like 
~ , being A'^ the volume of the substrate and 9 ranging from 1 — log 2/ log 3, for central target(s), 
to 1, for a single peripheral target. 



PACS numbers: 05.40.Fb,89.75.Hc 

I. INTRODUCTION 

The importance of first-passage phenomena stems from 
their fundamental role in stochastic processes that are 
prompted by a first-passage event; examples range from 
diffusion-limited growth, to neuron firing, and to the 
triggering of stock options [l]. These problems can be 
successfully dealt with by mapping the problem into a 
random-walk (RW) moving on a proper substrate and in 
the presence of an absorbing boundary, in such a way 
that first-passage properties can be expressed in terms of 
survival probability and mean hitting time. The relation 
between the first-passage probability and the probability 
distribution for a RW has been derived in different frame- 
works (continuous and discrete space and/or time), see 
e.g.[IHl]- In this context, the mean first-passage time 
(MFPT) to a target point, possibly averaged over all 
possible starting points, has been widely used in order 
to measure the transport efficiency jSl E] • 

The case of regular [2] and fractal lattices [7J [S] as 
substrate has been throughoutly investigated, and re- 
cently complex structures, especially (pseudo) scale-free 
networks, both random and deterministic, have also at- 
tracted a lot of interest [5HT^. The reason is that 
such networks are able to model some features typi- 
cal of many real systems and they can exhibit interest- 
ing, non-intuitive first-passage properties, whose under- 
standing could pave the way to engineered devices (see 
e.g. 1131 [H])- In particular, several examples were shown 
where the MFPT to the most connected node displays a 
sublinear dependence on the size N (9l . 

While previous works dealing with complex structures 
mainly focused on the case of one single trap placed on 
a central position, here we consider a special realization 
of deterministic scale-free structure where we study the 
MFPT for different arrangements of targets. Our re- 
sults, mainly rigorous, allow to evidence how deeply the 
position and/or the number of target (s) affect the first- 
passage properties. Indeed non-trivial effects are found 
especially when considering large and infinite graphs, 
which are used to describe macroscopic systems in the 



thermodynamic limit. More precisely, for the cases ana- 
lyzed we find that the the MFPT scales like ~ (up 
to logarithmic corrections), where N is the number of 
nodes and is a proper exponent. Moreover, ranges 
from 9o — 1 — log 2/ log 3 « 0.37 (when one single tar- 
get is placed on the central node) to 1 (when the target 
is placed on a peripheral node), which is expected to be 
the upper bound ^16j . Interestingly, the former sublin- 
ear scaling is recovered also when a number A^i°s2/iog3 
of targets is "centrally arranged". On the other hand, 
such an effective performance is by far not representative 
of the whole network and the MFPT averaged over the 
target site, denoted as f, is found to scale linearly with 
N. 

The paper is organized as follows: in Sec.|lT]and Sec. |III| 
we briefly resume the topological properties of the graph 
under study, while in Sec. |IV| we analyze the average 
time f. Then, we consider special arrangements of traps 
distinguishing between "central" , in Sec. [Vj and "periph- 
eral" , in Sec. |VI[ Finally, Sec. |VII| contains conclusions 
and discussions, while in Appendix [X]we have collected 
the details of analytical calculations. 



II. THE DETERMINISTIC SCALE-FREE 
GRAPH 

The graph under study can be built recursively in such 
a way that at the gth iteration we have the graph of 
generation g, denoted as Gg (see Fig. [Tjfor 5 = 3); here 
we briefly describe its main topological properties, for 
more details we refer to [TTl [T71 - fT^ . 

Starting from the so-called root constituted by one sin- 
gle node, at the first iteration one adds two more nodes 
and connects each of them to the root; the resulting chain 
of length three represents the graph of generation g = 1. 
We call Bi the set of sites introduced at the first genera- 
tion and linked to the root. Then, at the second iteration 
one adds two copies of Gi, namely G'l and C/"; for each 
of them one can distinguish the sets B[ and B'{, which 
make up the set B2, i.e. B2 = B'l U B'^; each node in B2 
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FIG. 1: (Color online) Deterministic scale-free graph of 
generation g = 3. According to the labeling used = 
{14, 15, 17, 18, 23, 24, 26, 27}. Moreover, the main hub has la- 
bel 1, while hubs of generation 2 have label 10 and 19; anal- 
ogously, nodes 11, 12,20,21 are rims of generation 1. 



is then directly connected to the root. Proceeding analo- 
gously, at the g-th iteration one introduces two replica of 
the existing graph of generation g — 1 and connects the 
root with all the sites belonging to Bg, henceforth called 
"rims" . 

This way, at the g-th generation the total number of 
nodes is Ng = 3^, the total number of rims is \Bg\ =2^ 
and the root turns out to be the main hub with coor- 
dination number 2(2^ — 1). Apart from the main hub, 
one can detect "minor hubs" , namely the two hubs per- 
taining to the subgraphs Gg-i and Gg-i] more generally, 
one can detect 2^ so-called hubs of generation g — 2 per- 
taining to the two copies of Gg-2 ^^"^ 
in cascade. Analogously, we will refer to rims of gener- 
ation g — n, with n = 0,1, ...,g — 1 (see Fig. [I]). Now, 
hubs follow a power-law degree distribution: P(k) ^ k^^ 
with exponent 7 — log 3/ log 2 « 1.59 [TT, while the re- 
maining nodes follow an exponential degree distribution: 
P{k) - (2/3)'^ = e-^^ where 7 = log(3/2) « 0.405 [18]. 

Finally, the average degree is {k)g = '}ZkP{k)/Ng = 
4[1 — (2/3)f], namely it approaches 4 as g — )■ 00; on the 
other hand, the second moment (fc^) is divergent [15] . 



III. VAN HOVE SPHERE AND AVERAGE 
DISTANCES 



main hub [T3]. Passing to the corresponding generat- 
ing functions, namely X{s,g) = J^'kLo -^(^^ d)^^ ^^'^ 
Y{s,g) = J2'^^oY{k,g)sS, we get 

Xis,g)^X{s,g-l) + 2Y{s,g-l)s (3) 
Y{s,g)=X{s,g~l)s + 2Y{s,g-l). (4) 

Now, by combining Eq. [3] and Eq. |4] we have 

l(s,5 + 2) =3X(s,.g + l)-2(l-s2)X(s,5). (5) 

The solution of this finite difference equation can be ob- 
tained by recalling the initial condition 



X{k,l) 



X(fc,0) = 4,o ^l(s,0) = l (6) 
4,0 + 4,1 ^l(s,l) = l + 2s, (7) 



due to the fact that the root of the graph (generation 
zero) is just the main hub, while the first generation 
graph is made up by the main hub and by two rims only. 
Hence we get [20] 
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which provide the cardinality of the van Hove surface of 
even and odd radius, respectively. Interestingly, these ex- 
pressions clarify the peculiar topology of the determin- 
istic scale free graph: due to the non monotonicity of 
X{k,g) and to the large number of surface points at a 
given distance, quantities obtained from the limit of a 
sequence of larger and larger graphs could present some 
anomalies [201. 

From X(s,g) we can also derive the average distance 
from the main hub dn- Let us define dij the average 
(chemical) distance from i to j, then we have 



N„ 



(10) 



For a graph of generation g, being X{k, g) the number 
of nodes at a distance k from the main hub and, analo- 
gously, Y{k, g) the number of nodes at a distance k from 
an arbitrary rim, the following equations hold: 

X[k,g) = X{k,g - 1) + 2Y{k -l,g- 1) (1) 
Y{k, g) = X{k -1,5-1) + 2Y{k, g-1). (2) 

Notice that X{k,g) represents the cardinality of the 
borders of the generalized Van Hove sphere with ra- 
dius k in the chemical distance and centered in the 



where H is the node representing the main hub. Now, 
the following relation holds 



(11) 



Then, from Eq. [5] we get 

—X{s,g + 2) = 3—Xis,g + l)+4sXis,g) 
OS OS 

- 2{1- s')^^X{s,g), (12) 
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which, for s = 1, gives 
Ng+2dH{9 + 'i)^iNg+^dH{9 + l)+AX{l,g) = Q. (13) 

From Eq. [5) X{l,g) = 39 TVo follows, as expected, and 
by plugging X{\,g) into Eq. Il3l we obtain the following 
recursive relation 



is taken over all possible couples. As shown in [TBI I23j . 

fg is subject to the following, rather strict, constraints: 



N„ 



<rg< 



Ng{k)dg 



iVglog 



N. 
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(20) 



where we used Eq. 18 for dg [H] 



g+2 



9d 
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whose solution is 
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(2g + l), 



(14) 



(15) 



where we used the initial condition di — 2/3. The aver- 
age distance to the main hub can also be expressed as a 
function of the volume, recalling that g = log A^g/log3 



log 3 



1 



log Ng 



(16) 



Such an average is characteristic of a network, in the 
sense that it provides an overall measure of the transport 
efficiency on that network |10j . but, on the other hand, it 
conceals the peculiarities of the graph itself and it does 
not allow to see whether there exists any target position 
which results to be particularly effective or ineffective. 

In order to deepen this point we start by considering 
the set {tij}i.j=i,...,N '■ For an arbitrary connected couple 
the mean time Uj can be calculated numerically by 
exploiting a method based on the so-called pseudoinverse 
Laplacian [rail^. 



While dg provides information about the "accessibil- 
ity" of the main hub itself, the average distance dg cal- 
culated over all pairs of nodes, provides a more global 
information, being tightly connected with dynamics phe- 
nomena on graphs such as spreading, random walks and 
synchronization [H [211 122] • The average distance is de- 
fined as 



V^" d 

''^=NgiNg-iy 



(17) 



The exact expression of dg for the graph under study has 
been obtained in |19| and is given by 



- 8iV,logjV, 
' ~ HNg 



l)log3 



^og Ng, 



(18) 




where the last expression holds in the thermodynamic 



limit N„ 



00. 



FIG. 2: (Color on line) Mean times tij divided by the volume 



The consistency between the leading behaviors of dg for graphs of generations 5 = 4, 5, 6, as indicated by the 



and dg, see Eqs. 



16 



and 



18 suggests that, despite its 



centrality, the main hub is not, in the average, spatially 
closer than any other site. In fact, the maximum distance 
among nodes grows logarithmically with the volume, that 



IS, maxij=i^...^Ar, 
rather compact. 



di 



2g, and the graph is generally 



IV. TOTAL MEAN FIRST PASSAGE TIME 

The average distance dg discussed in the previous sec- 
tion allows to determine an upper bound for the total 
mean first passage time fg defined as 



legend; only non degenerate data have been depicted. Notice 
that data are plotted in ascending order and versus a proper 
normalized interval. 



The intrinsic symmetry of the graph yields degener- 
ate values for tij which have been omitted in the plot of 
Fig. [2] From the plot one can notice that tij spans a wide 
range ^ 0{Ng) and the variance of the set {tij} (see 
the inset of Fig. IV I is shown to grow asymptotically as 
the square of the graph volume. The large variability of 



1 



NaiNg-l) 



(19) 



where tij is the mean time taken by a simple RW on Qg 
to first reach the site j starting from i, and the average 



tij mirrors the strong inhomogeneity of the graph under 
study; this is a further confirmation that the efficiency of 
first-passage processes is sensitively affected by the target 
position. 

From the numerical data collected we are able to cal- 
culate fg according to Eq. 



19 



as shown in Fig. IV (main 



figure) its asymptotic behavior is consistent with the lin 



ear law t„ ^ 3iV„, in agreement with Eq. 20 
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FIG. 3; (Color on line) Total mean first passage time fg (main 
figure) and related variance (inset) as a function of the 
graph volume Ng. Data points have been fitted with the 



power law < 



Ng and 



Ng , respectively. 



V. TARGET ON CENTRAL NODES 

As shown in [11], for a simple RW on Qg, the mean 
time to first reach the main hub starting from any of the 
rims in Bg is 



Ro 



(21) 



while the mean time to first reach any rim starting from 
the main hub is 



1, 



(22) 



where the mean is taken over all possible paths. Such 
results were used to calculate the mean time to first 
reach the main hub Tg{H) (where the mean is taken over 
all possible paths and over all possible starting points), 
which reads as 



1 



39-1 



4g)3S 



a 



-1/7 



32 /r 

(23) 

being 1 - I/7 = 6*0 = 1 - log 2/ log 3 « 0.37. Notice that 
the asymptotic behavior of Rg, Hg and Tg{H) is the same, 
namely (3/2)^. The exponent ^0 is very small if compared 
with the result found for fg and also with similar results 
found in the literature [SHTTl IT51 E5] , making the main 
hub a particularly effective node where to place a target. 

We now calculate the mean hitting time Tg{Bg) at any 
rim, averaged over all possible starting points ^ Bg and 
over all possible paths: 



1 



\B„ 



■E 



(24) 



By properly combining Hg and Rg we can write the fol- 
lowing recursive equation: 



rg{Bg) 



1 



,[3^-ii7, + (3^-i-l)r,_i 



2r,_i(6,_i)(3«-i-29-i) 



(25) 



In fact, if the RW starts from a site belonging to the 
subgraph Gg-i containing the main hub, then, in order 
to arrive to rims, it must pass through the main hub 
itself (in a mean time Tg^i{H)) from which rims are first 
reached in a time Hg. Conversely, if the starting point 
is in g'g_^ U ^g„i, the MFPT to rims is just Tg_i{Bg-i). 
The recursive relation can be solved to obtain 



39 

9(39 - 29) 



20 



18 



(26) 

Interestingly, the leading behaviour is the same as Tg(H), 

even if in this case the number of targets is 29 ~ Ng^'' . 
It is rather intuitive at this point that when targets are 
placed on the main hub and on rims the leading be- 
haviour of the hitting time is conserved. Indeed, one 
has 



TgiHUBg) - 



1 



39 - 29 -1 



[(3^-1 - l)r, 



r3_i(S3_i)(2.3f-i-29)]. (27) 



We finally consider the case of one single target placed 
on a rim, i.e. on a site b £ Bg] the calculation has been 
performed numerically exploiting the following expres- 



1 ''^ 

9 1=1 



i=l j = l 



where L is the Laplacian matrix of the underlying struc- 
ture. In the range considered, i.e. g € [0,6], data points 
can be fitted by the following 



Taib) 



N„ 



logWg' 



(28) 



which is consistent with the analytical result found in |16j 
(although the average here is taken over the fiat distribu- 
tion). The MFPT's Tg{H), Tg{Bg) and Tg{b) are shown 
and compared in Fig. [4] 



VI. TARGET ON PERIPHERAL NODES 

After having analyzed particularly effective arrange- 
ments for targets, i.e. main hub and rims, we now want 
to focus on trapping at the farthest nodes from the main 
hub which are expected to be rather difficult to reach. 

We notice that in a graph of generation g, the radius 
of the outmost van Hove surface centered on the main 
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this is found to be 



FIG. 4: (Color on line) Mean hitting times at the main hub 
Tg{H), at rims Tg{Bg) and at one rim Tg(b), as a function 
of the number of nodes Ng for a simple RW moving on Qg. 
Data from numerical calculations coincide with those calcu- 



lated analytically (for the former to cases see Eq. 23 Eq. 26 1 



and are compared with the leading behaviors ~ Ng 
~ Ng/ log Ng (continuous line). 



1/7 



and 



hub, i.e. the largest chemical distance achievable from 
the main hub, equals g and its cardinality is X{g, g) — g 
if g is even, otherwise X(g, g) — g + 1; such farthest nodes 
correspond to hubs {g even) or rims {g odd), of subgraphs 
of generation 1. We call Tg the set of nodes at a chemical 
distance g from the main hub so that \Tg\ = X{g,g). 

We will first consider the case of Tg working as a 
boundary, that is each node in Tg is occupied by a target, 
hence, exploiting the intrinsic symmetry of the arrange- 
ment, the problem will be solved analytically. Later we 
will focus on the case of one single target fixed on a node 
belonging to Tg and we will solve the problem numeri- 
cally. 



A. Multi-target 



3, 



(29) 



which exhibits the same asymptotic behavior of Hg and 
and [221 [11 . Then, the RW has to 



21 



of Rg (see Eqs. 
"bounce" from rim to hub and from hub to rims of inner 
and inner subgraphs, the pertaining MFPT's read as 



Rg_n = 4 X 2"/2 



and 



Hg^^ = 4 X 2"/2 



3X3 



3X9 
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(30) 



(31) 



We notice that i?g_„ and Hg^n (Eqs. |30|and|31[) display 



the same asymptotic behavior, namely (3/2)^2^^, which, 
interestingly, is the same behavior found for Tg{H) but 
"corrected" by a factor which depends on n as 2"/^: 



Rg-r 

Hg-n 



- 2"/2 Rg 

2"/^ H„. 



(32) 
(33) 



Some remarks are in order here. Although, for a given 
subgraph, the distance between the hub and a rim is 
the same and equal to 1, the mean time taken by the 
RW depends significantly on how "deep" the subgraph 
is (i.e. how small g — n is). Such a result is somehow 
counterintuitive because the probability to go directly 
from rim to hub [l/((? — n -\- 1)] and from hub to rim 
[2ff-"-i/(2S-» _ 1)] grows with n. Actually, what really 
matters is that if the right way is missed at the first step 
then it is more and more difficult for the RW to retrieve 
it as n gets larger. The factor 2"/^ is a clear indication of 
the fact that targets placed on inner nodes are in general 
rather difficult to reach. 

By merging the previous results one can find the aver- 
age time t^^n to reach a node in Jg, starting from the 
hub of any subgraph of generation g ^ n 



The basic step to calculate the hitting time Tg{Fg) 
when each node in Fg is occupied by a target, is to ob- 
tain Hg-n and Rg-n, which represent, respectively, the 
MFPT from the hub of a subgraph of generation g — nto 
any rim of the same subgraph, and the MFPT from any 
rim in a subgraph of generation g — n io the pertaining 
hub; for example, referring to Fig. [l] R^-i is the mean 
time to go from node 14 to 10, while iJ3_2 is the mean 
time to reach 11 starting from 10 (and equivalently for 
similar nodes). The detailed calculations are reported in 
the Appendix, while here we just summarize the main 
passages and results. 

We first consider the MFPT from the main hub to Fg] 
starting from the main hub the RW has first to reach the 



9-2 



rims in a time H„ (see Eq. 22 1 , from Bq the RW has to 



reach the hub of the subgraph of generation <? — 1 and 



(34) 



or from any rim of generation g — n 



9-2 9-2 

^f-"(-^9)=E^f-'+ E ^a-i- (35) 

/—I l—g—n—l 

These expressions are used in Appendix|A]to get the aver- 
age time Tg{Fg) to reach an arbitrary node f £ Fg, where 
the average is meant over all possible starting points ^ J-g 
and over all possible paths, namely 



Ng-lTgl^'-'- 
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Hence, we finally get 

^' 45(39 -2f) L V2/ 

- (45 X 29/2+2 _ 255 - 527) - 144 

^ ^^29/2^^;-^^ (37) 

where the last line contains the asymptotic behavior valid 
in the limit Ng — )■ oo (see Figj5]). Although the number 
of targets is ^ log Ng , the mean time to first reach any 
of them grows much faster than the mean time to reach 
the main hub Tg{H) ~ (3/2)9, g^iU slower than the 
average hitting time fg ^39. Indeed, recalling 1 — 5^ = 
1 - « 0.6845, we have that 

2 log 3 ' 

rgiJ'a) ^ fl'^ - Tg{H)'+^, (38) 
where 1 + 1/(27-2) w 1.85. 




FIG. 5: (Color on line) Mean hitting times at peripheral nodes 
Tg{J^g) (•) and Tg{f) (a), as a function of the number of nodes 
Ng for a simple RW moving on Qg . Data from numerical calcu- 
lations coincide with those calculated analytically from Eq.|26| 
and are compared with the leading behavior ~ Ng~^^^^'''' and 
~ Ng (continuous line), according to Eg. |38| and Eg. |39[ re- 
spectively. 



Results found for several generations are shown in Figj5] 
By means of fitting procedures we find that when the size 
is large, the mean hitting time Tg{f) scales like 

(39) 

where the exponent 0/ is consistent with the value 1, in 
agreement with [16j . 



VII. CONCLUSIONS 

In this paper we studied the mean first-passage prob- 
lem for a random walk on a deterministic scale- free struc- 
ture, focusing the attention on different arrangement of 
targets. First of all, we highlighted that, due to the high 
inhomogeneity of the structure, the mean first-passage 
time (MFPT) tij from a node i to j sensitively depends 
on the couple chosen, spanning a range ^ 0{Ng), 
being Ng the size of the structure. 

The case of a single target placed on the main hub 
turns out to be strategic as it gives rise to a slow sublinear 
scaling N^« , being 6*0 = 1 - I/7 = 1 - log 2/ log 3 w 
0.37. Interestingly, the same scaling is obtained in the 
case that the whole set of rims Bg, with cardinality ~ 

0(7Vg/^), is assumed as absorbing boundary. On the 
other hand, when only one single node b G Bg works as 
a target we get a faster growth ~ Ng/ log Ng, where the 
logarithmic correction accounts for the large coordination 
number of rims (see '1(5]). 

When peripheral positions for traps are considered, 
sublinear scaling are achievable provided that the num- 
ber of targets is large enough; in particular, when the 
absorbing boundary is given by the whole set of periph- 
eral nodes J^g with cardinality log Ng/ log 3 we have that 

the mean hitting time scale like ~ Ng~^^''^'^\ Conversely, 
when one single node f Bg works as a target our fits 
suggest a linear growth ^ Ng. 

Therefore, the transport efficiency sensitively depends 
on the number and on the position of targets, moreover, 
the existence of very "close" couples, such as main hub 
and rims, is evened off by rather "distant" couples, giving 
rise to a mean hitting time fg, averaged over all possible 
pairs, which scales linearly with Ng. 



Appendix A: Mean FPT on Tg 



B. Single Target 

We now turn to the case of one single target placed on 
a node f G Tg, which has been investigated numerically 
according to 

N N N 

^.(/) = ^Et./ = ^EE(-L-).. 

^ 2—1 ^ i — 1 7 — 1 



In this appendix we report the detailed calculation of 
the mean time taken by a simple RW on a graph Qg 
and started from the hub of a subgraph of generation 
n, to first reach any rim in the same subgraph, which 
is denoted as and the MFPT from any rim in a 
subgraph of generation n to the pertaining hub, denoted 
as Rn- Such times will then be merged in order to obtain 
the mean time Tg{J-g) (averaged over all possible starting 
points and over all possible paths) to first reach any node 
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belonging to the set Fg and meant to be occupied by a 
target. 

Let us start from the highest generation, i.e. n = g. 
As for Hg, this has aheady been calculated in [TT] (see 
Eq. 22 ) . Now, from rims of generation g, in order to reach 



J-"g, we first need to pass through the hub of gncration 



5- 



1 and we can write 



Rg-i - 



9-2 



(tf + l + Rg) + l + j2a + tf 



Rn 



.(Al) 



1=1 



In fact, starting from any rim of generation g, a RW has 
g options: it goes to the main hub and then returns to 
rims in a time iJg, either it goes directly to the hub of 
generation .g — 1, either it reaches a hub of generation 
g — I from which it return to rims in an average time Hi . 
By plugging Eqs. 22 and 21 in Eq. ( Al ) we get: 



Rq-l 



28 

y 



(A2) 



Now, from the hub of generation g ~ 1 the RW needs to 
reach the pertaining rims and this is first accomplished 
in a mean time given by: 



H. 



1 



9-1 



2(29-1 _ 1) [ 

(g-i)-2 

{l + Rl+ Hg^l) + (1 + + 



(A3) 



in fact, the RW can jump directly in any of the 
2(2^9-1)-!) rims of generation {g — 1) — 1, either it can 
reach rims of g — 1 from which it needs to return to the 
hub of generation {g — 1) in a time Rg-i (Eq. Al), either 
it can jump to a rim of the subgraph of generation g — I 
from which it can return to the hub in i?/. 

As can be deduced from the previous equations, Hg_n 
and Rg-n are deeply connected with each other, such an 
interplay allows to build up a system of recursive equa- 
tions. We can in fact write the following general equa- 
tions 



9 « 2(25-" - 1) 



2(9-")- 



(A4) 



(g-n-2) 

i^ + Rl+ Hg-n) + (1 + Rg-n + Hg^, 

1 = 1 



and 



Rg—n 



k-ffg-n + 1 + Rg-n) 

g - nl 

9-2 

J2{^+Hl+Rg-r 
1 = 1 



(A5) 



Now, by plugging in Eqs. "KE and "KA the expressions 



Eqs. 21 and [22} with some algebra we get 

'gX 9-" 

Rg-n = Hg-n+1 + 2 -|- 4 ' ' 
Hg-n — 2i?g_„+i + 1 



2, 

3X9-n 



9 V2 



, (A6) 
(A7) 



By properly combining these equations we have 

28 /SX^ " 

^9-n = ^9-"+2-l+y ( 2) (A8) 



and using Eq. 21 we have 



(A9) 



Repeating the same calculations for Hg^n we obtain sim- 
ilar results, namely: 

i?,_„ = 4 2-/2 - (2^ - 3) - ^ Qy^" . (AlO) 

Following analogous operations for n odd, we reach 
the same result provided that n is replaced by n 4- l.We 
notice that, for a given generation g, the calculation of 
Tg{J-) involves either Rg-2n or Rg-2n+i occur according 
to whether g is even or odd. Hence, in the following, 
without loss of generality, we focus on the case of t; — n. 

From Hg-n and Rg-n it is now possible to derive a 
closed form expression for the trapping time at J^g. We 
now use the partial results found so far to finally calculate 
Tg{J^g). First of all, let us calculate the first passage 
time (J-g) on any site in for a RW started from the 
main hub; when g is even, J-g is made up of nodes which 
correspond to hubs of generation 1, vice versa, when g is 
odd Tg is made up of nodes which correspond to rims of 
generation 1. Let us now focus on g even; we then reach 
the following expression: 



9-2 



tfi^g) - Y ^3-1 + H, 



g-l, 



(All) 



which can be straightforwardly generalized to the case of 
a RW started on a hub of generation g — n 

s-2 

C(-^9)-E(^^'-'+^9-')' (A12) 
I 

from which we get 

tf_„(^g) = 8(2^/^-2-/2) 



32 

~5 



g-n 



- 1 



2(g~7i).(Al3) 



8 



Let us now calculate the mean time t^{Fg) taken by a 
RW started from a rim of generation g — n io first reach 
any node in Tg] the following relation holds 



9-2 9-2 

l—l l—g—n—1 



(A14) 



which is just a truncated version of sum (A12), so that 
we get 



1 



56 /3 



15 V 2 



+ 2{g -n) + 



17 



(A15) 



In order to simplify the notation we introduce T„ which 
represents the mean time taken by the RW to first reach 
a site / in J-g and started from a node at a chemical dis- 
tance n from /. Since we assumed g to be even, when 
n is even (odd) the starting point is a hub (rim). By 
exploiting the intrinsic symmetry of the graph, T„ can 
be written as the sum of two contributes, according to 
whether the starting node belongs to the subgraph Gg-i 
(and therefore the RW has to pass through the main hub 
before reaching the target) or to G'g-i U Gg^i (and there- 
fore the RW has to pass through Bg); for subgraphs of 
smaller generation it holds analogously. The first con- 



tribute reads as 



9-1 



J2 2^(3^-" - l)(rf_„ + T,_„+i) + T,_„+i, (A16) 



n=l 
n odd 



in fact the RW has first to reach the main hub of the 
pertaining subgraph in a mean time Tg_„, from which it 

reaches J-g in a time Tg^n+i- The factor 2 2 accounts 
for the number of subgraphs of generation g — n present. 
The second contribute reads as 

9-2 

4 J2 2^(3f-"-2^-")(rji„+T,_„+i) + 2^""7^9-n+i. 

n even 

(A17) 

where we considered RWs starting from nodes such that 
in order to reach Tg has first to reach the rims of the 
pertaining subgraphs or Bg itself. We sum together the 
previous expressions to obtain 



^9 (-^9 



29^ 



45(39 - 23) 

3X9 



^'2 



29 



(45 X 2^/2 - 97) 



(45 X 29/2+2 „ 25.9 - 527) - 144 



for g odd analogous results are expected. 
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